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Abstract 


Command  procedures  have  been  developed  for  the  Army  Research 
Laboratory  molecular  beam  epitaxy  (MBE)  computer  control  system 
that  allow  a  user  to  automatically  outgas  and  desorb  the  oxide  from 
substrates  before  growth,  as  well  as  set  substrate  temperatures 
based  on  pyrometer  readings  during  growths.  These  procedures  allow 
completely  unattended  growth  of  structures  once  suitable 
temperatures  have  been  determined. 
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Figure 


i.  Temperature  as  a  function  of  time  during  growth  of  two  consecutive 
PHEMT  wafers  .  4 


Background 


The  Army  Research  Laboratory  molecular  beam  epitaxy  (MBE)  system  is 
controlled  by  a  PC -based  system  that  uses  the  "Molly"  software  package 
supplied  by  EPI.  The  Molly  system  provides  a  script  language  that  can  be 
used  to  create  command  procedures  that  execute  customized  sequences  of 
actions  on  the  MBE  system.  Allowed  actions  are  reading  and  setting  cell 
temperatures,  opening  and  closing  shutters,  reading  pressure  gauges,  set¬ 
ting  the  azimuthal  rotation  speed  of  the  substrate,  and  turning  the  substrate 
holder  to  the  growth  and  load  positions. 

Before  deposition  of  material  is  started,  the  native  oxide  on  the  substrate 
must  be  desorbed.  This  is  done  by  heating  the  substrate  to  a  predetermined 
temperature  under  As  overpressure.  Due  to  differences  in  emissivity  of  the 
various  substrate  holders  because  of  design,  size,  and  age,  a  specific  ther¬ 
mocouple  setting  cannot  be  expected  to  produce  the  desired  temperature 
with  enough  precision.  Hence,  small  corrections  to  the  set  points  must  be 
applied,  based  on  reading  of  the  optical  pyrometer. 

Many  device  structures  consisting  of  various  material  combinations  may 
require  different  substrate  temperatures  during  the  growth  sequence.  For 
the  reasons  mentioned  above,  it  is  not  possible  to  preprogram  a  tempera¬ 
ture  profile  for  the  thermocouple  since  it  may  not  result  in  the  desired  tem¬ 
peratures  as  read  by  the  pyrometer.  An  operator  must  therefore  be  present 
to  make  the  final  temperature  adjustments.  Automation  of  these  adjustments 
leads  to  time  savings  and,  in  some  cases,  greater  precision  and  reproduc¬ 
ibility. 


Script  for  Outgassing 


The  logic  behind  the  script  developed  for  substrate  outgassing  is  as  fol¬ 
lows:  The  substrate  temperature  is  first  raised  to  a  thermocouple  set  point, 
at  and  above  which  the  pyrometer  operates  (440  °C).  Below  this  tempera¬ 
ture  the  bandgap  of  GaAs  is  large  enough  that  the  substrate  is  transparent, 
which  results  in  a  very  high  pyrometer  reading  since  the  pyrometer  "sees" 
the  heater  elements  behind  the  wafer.  After  the  minimum  operation  tem¬ 
perature  has  been  reached,  the  computer  starts  reading  values  from  the  py¬ 
rometer.  (The  critical  temperature  is  considered  reached  when  the  pyrom¬ 
eter  reading  exceeds  450  °C  and  the  time  derivative  of  the  temperature  is 
positive.)  The  pyrometer  can  be  read  from  either  the  analog  or  digital  out¬ 
puts  of  the  instrument.  The  digital  output  is  greatly  preferred  due  to  the 
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much  lower  noise  level  on  the  signal.  (The  analog  output  was  used  origi¬ 
nally  because  of  a  malfunction  with  the  digital  output  port  and  produced  a 
signal  that  required  heavy  averaging  and  thus  more  script  code.) 

The  ramp  between  the  minimum  operating  temperature  of  the  pyrometer 
and  the  final  outgassing  temperature  is  approximated  with  a  10-step  stair¬ 
case  during  which  a  PID-control  (proportional  integration  and  differential) 
routine  drives  the  pyrometer  temperature  to  each  set  point.  When  the  final 
outgassing  temperature  is  reached,  it  is  held  constant  with  software  PID 
control  for  a  specific  time,  typically  10  min.  After  this,  the  temperature  is 
ramped  down  to  a  desired  growth  temperature  with  the  same  control  algo¬ 
rithm.  The  outgassing  script  has  now  finished  and  a  growth  script  can  be 
started.  The  code  for  the  outgassing  sequence  is  called  Genoutg.cmd  and 
can  be  used  for  both  GaAs  and  InP  substrates.  A  parameter  flag  must  be  set 
to  provide  information  about  which  substrate  is  used.  The  control  code  for 
the  outgassing  script  is  shown  in  in  the  appendix  (Command  File  Listings). 


Script  for  Temperature  Control  During  Growth 


By  running  the  script  Tsub.cmd  in  parallel  with  the  structure  recipe  during 
growth,  one  can  change  the  substrate  temperature,  as  read  by  the  pyrom¬ 
eter,  from  the  structure  recipe.  Again,  PID  software  control  is  used.  The  tem¬ 
perature  set  point  is  stored  in  a  global  variable  (see  next  section),  which  is 
read  every  5  s  by  Tsub.cmd  and  can  be  changed  in  the  recipe  by  loading  of 
command  files  named  TsubXXX.cmd,  where  XXX  is  the  temperature.  Al¬ 
lowed  values  are  450  °C  <  XXX  <  700  °C  in  5°  increments. 


Starting  the  Scripts 


The  appendix  shows  an  example  of  a  structure  recipe  for  a  PHEMT 
(pseudomophic  high-electron-mobility  transistor),  which  will  be  used  to 
illustrate  the  start  of  the  scripts.  All  command  procedures  are  invoked  with 
the  Molly  command  load.  A  few  comments  about  this  command  are  appro¬ 
priate  at  this  point.  The  Molly  language  also  includes  a  command  nsh,  which 
can  "launch"  another  command  procedure.  The  difference  between  load¬ 
ing  and  launching  is  that  a  loaded  command  file  must  execute  to  comple¬ 
tion  before  the  next  step  in  the  calling  (loading)  file  can  continue.  A  launch, 
on  the  other  hand,  starts  the  new  command  file  and  immediately  continues 
with  the  next  step.  In  other  words,  the  calling  and  the  called  procedure  will 
run  in  parallel.  For  unknown  reasons,  the  Molly  spreadsheet  accepts  only 
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the  load  command  and  not  nsh.  Incidentally,  the  new  EPICAD  program  al¬ 
lows  both  load  and  nsh. 

The  first  command  file,  Stdlogl.cmd  (see  appendix),  starts  a  data  logger. 
This  is  a  personal  preference  of  mine  and  is  not  needed  for  the  pyrometer 
control.  The  same  is  true  for  the  next  command,  Rot5.cmd,  which  sets  the 
rotation  speed.  Next,  the  outgassing  procedure,  Genoutg.cmd,  follows.  This 
should  run  to  completion  before  the  growth  begins  and  is  therefore  started 
directly  with  a  load. 

Since  we  clearly  want  to  allow  the  temperature  controller  and  the  following 
structure  recipe  file  to  run  in  parallel,  we  must  load  an  additional  interme¬ 
diate  file,  Pyroc.cmd  (see  appendix),  which  contains  only  the  nsh  command 
that  launches  the  desired  file  (Tsub.cmd)  (see  appendix).  This  is  done  im¬ 
mediately  after  the  outgassing. 

The  temperature  set  point,  which  is  used  in  the  PID  control,  is  the  same  as 
the  last  temperature  in  Genoutg.cmd.  This  temperature  is  stored  in  the  glo¬ 
bal  variable  Tsubjiezv,  which  is  declared  in  the  file  C:\molly\lib\local.cmd, 
a  command  file  that  is  executed  automatically  when  the  Molly  system  is 
started.  One  can  make  subsequent  changes  to  the  substrate  temperature  by 
loading  files  called  TsubXXX.cmd,  where  XXX  is  the  temperature  (see  ap¬ 
pendix).  Allowed  values  are  450  °C  <  XXX  <  700  °C  in  5°  increments.  These 
command  files  only  set  the  global  variable  to  a  new  value,  so  their  execu¬ 
tion  time  is  negligible.  If  a  different  temperature  is  desired,  any  of  the  exist¬ 
ing  files  can  be  copied  and  renamed  to  the  new  temperature  and  the  set 
point  edited.  (Temperatures  below  450  °C  should  not  be  used  for  the  obvi¬ 
ous  reason  that  the  pyrometer  does  not  work  in  that  range.)  In  the  PHEMT 
recipe,  loading  Tsub500.cmd  lowers  the  temperature  from  600  °C,  which  is 
the  temperature  at  the  finish  of  the  outgassing,  to  500  °C  just  before  the 
InGaAs  layer  is  started. 

Before  the  InGaAs  growth  starts,  another  GaAs  layer  is  inserted  after  the 
Tsub500.cmd.  This  illustrates  an  important  limitation  of  the  PID  routine.  If 
the  temperature  change  is  too  large,  an  abrupt  change  in  set  point  may  re¬ 
sult  in  a  substantial  under-  or  overshoot.  An  undershoot  can  be  particularly 
problematic  since  it  may  drive  the  substrate  temperature  below  the  operat¬ 
ing  range  of  the  pyrometer.  To  avoid  this,  the  Tsub.cmd  file  always  applies 
a  ramp  with  a  rate  of  40° /min.  Hence,  in  the  PHEMT  example,  the  new  set 
point,  which  is  100  °C  below  the  previous  one,  will  not  be  reached  until 
150  s  after  Tsub500.cmd  is  loaded,  requiring  the  extra  GaAs  layer  (180  s  is 
used  in  the  recipe  for  extra  margin). 

When  the  recipe  has  finished,  the  Tsub.cmd  should  be  terminated.  This  is 
done  by  loading  Kill_pyr.cmd  (appendix).  The  latter  executes  the  command 
kill(pyro_id ,  SIGKILL)  in  which  pyro_id  is  another  global  variable  defined  in 
C:\molly\lib\local.cmd.  Similarly,  the  logging  is  stopped  with  the  file 
Kill_log.cmd  (appendix).  The  files  Rot0.cmd  and  Posl.cmd  turn  off  the 
azimuthal  rotation  and  turn  the  substrate  to  the  loading  position.  (Note 
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that  the  pyrometer  control  definitely  must  be  turned  off  before  the  substrate 
is  turned  in  another  direction.) 

Figure  1  shows  the  results  of  two  consecutive  growths  (run  numbers  1281 
and  1282)  with  the  automated  outgassing  and  the  change  of  substrate  tem¬ 
perature  during  deposition.  Despite  the  fact  that  two  different  types  of  sub¬ 
strate  holders  (EPI  UNI-Block  and  Varian  nonbonded)  were  used,  the  py¬ 
rometer  traces  are  virtually  identical  beyond  the  outgassing  phase. 


Figure  1.  Temperature 
as  a  function  of  time 
during  growth  of  two 
consecutive  PHEMT 
wafers.  Solid  trace 
represents  thermo¬ 
couple  data  and 
dotted  trace  represents 
pyrometer 
temperature. 


Sample  A1281 


Time  (s) 

Sample  A1282 


0  1000  2000  3000  4000  5000 
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Final  Comments 


The  files  Genoutgas.cmd,  Tsub.cmd,  TsubXXX.cmd,  and  Pyroc.cmd  are  all 
stored  in  C:\u\mbe\.  This  allows  the  use  of  the  file  name  only  (without  the 
explicit  path  name)  in  the  growth  recipe.  As  mentioned  above,  the  global 
parameters  are  created  in  C:\molly\lib\local.cmd. 

The  parameters  that  are  coded  into  the  script  files  all  have  comments  next 
to  their  declarations.  It  should  be  fairly  obvious  from  the  script  listings  what 
the  parameters  represent.  Some  optimization  has  been  done  and  the  present 
values  do  work.  However,  further  optimization  can  probably  be  done  and 
certainly  refinements  are  possible.  In  particular,  I  am  planning  to  allow  a 
step  change  in  temperature  if  the  difference  between  the  old  and  new  set 
points  is  below  a  certain  value.  Once  I  am  satisfied  that  the  outgassing  rou¬ 
tine  works  perfectly  for  both  InP  and  GaAs,  I  will  probably  create  separate 
files  so  that  the  operator  does  not  have  to  edit  the  parameter  SubToday. 


Acknowledgment 


I  used  the  PID  controller,  coded  and  optimized  by  Rich  Leavitt. 
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Appendix.  Command  File  Listings 


PHEMT8.wbl 

Genoutgas.cmd 

Stdlogl.cmd 

Pyroc.cmd 

Tsub.cmd 

Tsub550.cmd 

Kill_pyr.cmd 

Kill_log.cmd 


Molly  spreadsheet  recipe  for  a  PHEMT 
General  outgassing  procedure  for  GaAs  and  InP 
Starts  the  standard  logging  procedure 
Launches  Tsub.cmd 
General  PID  controller  for  MBE  growth 
Changes  the  pyrometer  set  point  to  550  °C 
Terminates  Tsub.cmd 
Terminates  the  data  logging 
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phemt8  | EDIT  FRAMED  PARAMETER | 


Appendix 


PHEMT8.wbl 
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Genoutgas.cmd 


Appendix 


/**/ 


/***************************************************************************************/ 
/*  */ 


/*  This  command  file  desorbs  the  oxide  from  an  substrate  which  can  be  either  */ 
/*  GaAs  or  InP  (not  tested  yet)  */ 
/*  Temperature  control  via  pyrometer  and  PID  routine  from  Rich  Leavitt  */ 
/*  */ 


/*  Stefan  Svensson,  ARL  Mar  9,  2000  */ 

/***************************************************************************************/ 


/**/ 


#include 

#include 

#include 

#include 

#include 

#include 

#include 

#include 

/**/ 

/*=  =  =  =  =  =  : 


<stdlib . h> 
<signal . h> 
<unistd. h> 
<cells . h> 
<gauges . h> 
<shutters . h> 
<time . h> 

<mbe . h> 


*/ 


/* 


Declarations  and  initial  values 


*/ 


/*====■ 

/**/ 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/**/ 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 


*/ 


/- 

- \  -TPoutgas 

*/ 

*/ 

/ 

\ 

*/ 

/ 

\ 

*/ 

TC_work- 

/-/ 

\ -  TPGrow 

*/ 

/ 

*/ 

/ 

*/ 

/ 

*/ 

/ 

*/ 

— / 

*/ 

*/ 

*/ 

*/ 

*/ 

Ramp_t ime 1 

Ramp_ 

_time2  Ramp_time3 

Wait_ 

timel  Outg_time 

|  Pyro  reading  starts 

*/ 

Substrate; 

/* 

Substrate  flag 

*/ 

GaAs 

= 

1; 

/* 

GaAs 

*/ 

InP 

= 

0; 

/* 

InP 

*/ 

SubToday 

= 

1; 

/* 

Today' s  substrate 

*/ 

TPoutgasInP 

= 

550.0; 

:  /* 

Pyrometer  Max  outgassing  temp  InP 

*/ 

TPoutgasGaAs 

= 

630.0; 

:  /* 

Pyrometer  Max  outgassing  temp  GaAs 

*/ 

TPoutgas; 

/* 

Pyrometer  Max  outgassing  temp  General 

*/ 

TP_growInP 

= 

525.0; 

:  /* 

Pyrometer  Growth  temperature  InP 

*/ 

TP_growGaAs 

= 

600.0; 

:  /* 

Pyrometer  Growth  temperature  GaAs 

*/ 

TP_grow; 

/* 

Pyrometer  Growth  temperature  General 

*/ 

TPdelta; 

/* 

Pyrometer  Difference  between  last  two  readings 

*/ 

TPlast ; 

/* 

Pyrometer  Previous  reading 

*/ 

TP sub; 

/* 

Pyrometer  general  substrate  temperature 

*/ 

TPO ; 

/* 

Pyrometer  reading  at  start  of  PID  control  ramp 

*/ 

TPwork 

= 

440.0; 

:  /* 

Pyrometer  min  real  temperature 

*/ 

TP step; 

/* 

Pyrometer  step  length  in  PID  control  ramp 

*/ 

TCsub; 

/* 

Thermocouple  general  substrate  temperature 

*/ 

TC_work; 

/* 

Thermocouple  Min  operating  temp  (first  ramp  target) 

*/ 

TC_workGaAs 

= 

650.0; 

:  /* 

Thermocouple  Min  operating  temp  GaAs 

*/ 

TC_workInP 

= 

500.0; 

:  /* 

Thermocouple  Min  operating  temp  InP 

*/ 

TC_test 

= 

450.0; 

:  /* 

Thermocouple  start  testing  for  pyrometer  operation 

*/ 

Ramp_rate; 

/* 

Thermocouple  Ramp  rate 

*/ 

Ramp_t ime 1 

= 

O 

o 

1 — 1 

/* 

Minutes  ramp  time  to  temp  where  pyrometer  works 

*/ 

Ramp_t ime 2 

= 

5.0; 

:  /* 

Minutes  ramp  time  to  outgassing  temp 

*/ 

Ramp_t ime  3 

= 

3.0; 

:  /* 

Minutes  ramp  time  to  growth  temp 

*/ 

Wait_t imel ; 

/* 

Minutes  to  wait  for  pyrometer  to  start  working 

*/ 

Wait_t ime2 

= 

0.1; 

:  /* 

Minutes  to  wait  between  pyrometer  readings 

*/ 

Wait_t ime3 

= 

0.5; 

:  /* 

Minutes  extra  wait  before  pyrometer  control 

*/ 
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Genoutgas.cmd 


Appendix 


double 

Outg_time 

=  10.0; 

/* 

Minutes 

outgas  time 

double 

Time_step; 

/* 

Time  step 

length  in  ramp 

double 

Numb_steps 

=  10.0; 

/* 

Number  of 

steps  that  ramp  is 

broken  into 

int 

i; 

/* 

Counter  for  pyrometer  readings 

int 

I_pyro 

=  5; 

/* 

Number  of 

pyrometer  readings 

to  average 

double 

to, 

control_time, 
Tstep  =  1, 
sigma, 
delta. 

Prop  =  72, 

Integ  =  4.7, 

Deriv  =  1.35, 

Tmin  =  25, 

Tmax  =  850, 
errtot=0 , 

TCsetp, 
oldT, 
newT, 
curr_err, 
prev_err ; 

string  sznewT, szct; 


/**/ 

/*===========================================================: 

/ k =============  initialize  parameters  for  the  right  substrate 

/*===========================================================: 

/**/ 

if(SubToday  ==  GaAs) 

{ 

TPoutgas  =  TPoutgasGaAs ; 

TP_grow  =  TP_growGaAs; 

TC_work  =  TC_workGaAs; 

} 

if(SubToday  ==  InP) 

{ 

TPoutgas  =  TPoutgasInP; 

TP_grow  =  TP_growInP; 

TC_work  =  TC_workInP; 


*/ 

*/ 

*/ 

*/ 

*/ 


*/ 

*/ 

*/ 


echo("  TPoutgas  =  ", TPoutgas); 
echo("  TP_grow  =  ",TP_grow); 
echo("  TC_work  =  ",TC_work); 

/**/ 

/*=========================================================; 

/ 'k =======================  Ramp  substrate  to  the  range  where 

/ k =======================  the  pyrometer  starts  working 

/*=========================================================; 

/**/ 


*/ 

*/ 

*/ 

*/ 


/*  Read  current  TC  set  point 


while (  TCsub  <=  0.0  )  TCsub  =  temp (  subs  ) ; 

Ramp_rate  =  fabs (  (TC_work  -  TCsub  ) /Ramp_timel ) ; 
set_ramp (  subs  ,Ramp_rate  );  /* 

set_temp (  subs  ,TC_work);  /* 

echo ("") ; 

echo("  Ramping  substrate  to  pyrometer  range  (Thermo-couple= 
echo("  This  will  take  " , Ramp_t ime 1 , "  min."  ); 

echo("");  /*  Wait  until  ramp  is  done 

echo("  Waiting  for  substrate  to  thermalize  and  pyrometer  to  start  working  ") ; 
Wait_timel  =  Ramp_t ime 1 * 1 . 1 ; 

sleep  (  Wait_timel * 60 ) ;  /* 


Ramp  thermocouple 

to  temp  where  pyro  works 


',TC_work,"  )"); 


/* 


Start  checking  pyro  just  before*/ 
ramp  finsihes  */ 


10 


Genoutgas.cmd 


Appendix 


/ 'k ===========  start  testing  for  pyrometer  operation  ================================ 

/*=================================================================================: 

/**/ 

echo("  Looking  for  pyrometer  stability  ") ; 

TPdelta  =  -1.0;  /*  initialize 

TPlast  =  1000.0;  /*  Search  for  positive  TPdelta 

TCsub  =  0.0; 

TPsub  =  0.0; 

while  (TPdelta  <0.0  | |  TPsub  <  TPwork) 

{ 

i=0 ; 

TPsub  =  0.0; 
while  (i  <  I_pyro) 

{ 

i=i+l ; 

TPsub  =  TPsub  +  reading (pyrometer ) ; 

/*  echo("  i,  TPsub  =",i, TPsub);  */ 

sleep (0.5) ; 

} 

TPsub  =  TPsub/I_pyro; 

TPdelta  =  TPsub  -  TPlast; 

TPlast  =  TPsub; 

/*  echo("  TPsub,  TPdelta  =  ",  TPsub,  TPdelta  );  */ 

sleep  (  Wait_time2*60)  ; 


*/ 

*/ 

*/ 


*/ 

*/ 


sleep  (  Wait_time3*60) ;  /*  Wait  an  extra  30  sec  (chicken!) 

/**/ 

/*==================================================================================== 

/ *=======================  Ramp  to  out gas  temp  ======================================= 

/*==================================================================================== 

/**/ 


echo("  Starting  pyrometer  control  and  upramp") ; 
echo ("") ; 

Time_step  =  Ramp_time2*60 . 0/Numb_steps; 
control_time  =  Time_step; 

TP0  =  reading (pyrometer ) ; 

TPstep  =  (TPoutgas  -  TP0 ) /Numb_steps; 

echo("  Current  temperature  =  ",TP0); 

/*  echo ( "control_time  =  ", control_time) ;  */ 

/*  echo ("TPstep  =  ", TPstep  );  */ 

/*  echo("Tstep  =  ",Tstep  );  */ 

i  =  0; 

while  (i  <  Numb_steps) 

{ 

i  =  i+1 ; 

newT  =  TPstep*i  +  TP0; 

echo ("New  temp : newT)  ; 

sigma  =  -setp (subs) *Integ*Prop/100; 

prev_err  =  reading (pyrometer ) -newT; 

set_ramp (subs,  1000)  ; 

oldT  =  newT; 


/*  Time  per  step  in  sec 


tm  =  t  ; 

tO  =  t; 

while  (  t-tO  <  control_time  ) 

{ 

sleep (Tstep  -  (t  -  tm) ) ; 

curr_err  =  reading (pyrometer ) -newT; 

sigma  =  sigma+Tstep*curr_err ; 

errtot  =  errtot+Tstep*curr_err*curr_err; 

delta  =  ( curr_err-prev_err ) /Tstep; 

TCsetp  =  - (curr_err+sigma/Integ+delta*Deriv) *100/Prop; 

if  (TCsetp  <  Tmin)  TCsetp  =  Tmin; 

if  (TCsetp  >  Tmax)  TCsetp  =  Tmax; 

set_temp (subs, TCsetp) ; 

prev_err  =  curr_err; 


*/ 


*/ 

*/ 

*/ 


*/ 
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Genoutgas.cmd 


tm  =  tm  +  Tstep; 
tm=t ; 


/**/ 

/*================================================================== 

/  =======================  Outgas  ================================== 

/ ^ ================================================================== 

/**  / 

echo("  Outgas  temperature  reached  ")  ; 
echo("  Outgassing  for  ",  Outg_time,  "  min"); 

/**/ 

control_time  =  Outg_time*60 . 0; 
newT  =  TPoutgas; 

/*  echo ("New  temp : newT) ;  */ 

sigma  =  -setp (subs) *Integ*Prop/100; 

prev_err=reading (pyrometer) -newT; 

set_ramp (subs, 1000) ; 

oldT  =  newT; 

tm  =  t; 

tO  =  t; 

while  (  t-tO  <  control_time  ) 

{ 

sleep (Tstep  -  (t  -  tm) ) ; 

curr_err  =  reading (pyrometer ) -newT; 

sigma  =  sigma+Tstep*curr_err ; 

errtot  =  errtot+Tstep*curr_err*curr_err; 

delta  =  ( curr_err-prev_err ) /Tstep; 

TCsetp  =  - (curr_err+sigma/Integ+delta*Deriv) *100/Prop; 

if  (TCsetp  <  Tmin)  TCsetp  =  Tmin; 

if  (TCsetp  >  Tmax)  TCsetp  =  Tmax; 

set_temp (subs,  TCsetp)  ; 

prev_err  =  curr_err; 

tm  =  tm  +  Tstep; 

tm  =  t; 

} 

/**/ 

/*================================================================== 

/ k =======================  Ramp  to  growth  temp  ===================== 

/ * ================================================================== 

/**/ 

echo  ("") ; 

echo("  Done  outgassing  ") ; 

echo("  Starting  " , Ramp_t ime 3 , "  min  down-ramp"); 
echo ("") ; 

control_time  =  Ramp_time3*60; 

TPO  =  reading (pyrometer ) ; 

echo ("Current  temperature  =  ",TP0); 


newT  =  TP_grow; 

echo ("New  temp :", newT) ; 

sigma  =  -setp (subs) *Integ*Prop/100; 

prev_err  =  reading (pyrometer ) -newT; 

set_ramp (subs,  1000)  ; 

oldT  =  newT; 

tm  =  t  ; 

tO  =  t; 

while  (  t-tO  <  control_time  ) 

{ 

sleep (Tstep  -  (t  -  tm) ) ; 

curr_err  =  reading (pyrometer ) -newT; 

sigma  =  sigma+Tstep*curr_err ; 

errtot  =  errtot+Tstep*curr_err*curr_err; 

delta  =  ( curr_err-prev_err ) /Tstep; 

TCsetp  =  - (curr_err+sigma/Integ+delta*Deriv) *100/Prop; 
if  (TCsetp  <  Tmin)  TCsetp  =  Tmin; 
if  (TCsetp  >  Tmax)  TCsetp  =  Tmax; 


:  =  =  =*/ 
■===*/ 
■===*/ 


*/ 

*/ 

*/ 
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set_temp (subs, TCsetp) ; 
prev_err  =  curr_err; 
tm  =  tm  +  Tstep; 

tm  =  t; 


/**/ 

/*======================================================================================*/ 

/*=======================Ready  to  grow  =================================================*/ 

/*======================================================================================*/ 

/**/ 

Tsub_new  =  TP_grow; 

echo("  Substrate  outgassed  and  oxide  desorbed  ") ; 
echo("  Ready  for  growth"); 
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Stdlogl.cmd 


/************************************************************************/ 


/*  */ 

/*  This  command  procedure  starts  a  data  log  of  temperatures,  */ 

/*  shutter  positions  and  the  flux  gauge  during  a  growth.  */ 

/*  The  data  is  stored  in  the  file  */ 

/*  C : /u/mbe/recipes/stef an/log_data/outf ile . dat  */ 

/*  */ 


/************************************************************************/ 

log_id  =  logger (20.0,  /*  log  every  20  seconds  */ 

't'  , 

'temp ( subs )  '  , 

' t  emp ( Ga ) ' , 

'temp  (A13)  '  , 

'temp  (A14 )  '  , 

'temp  (In)  '  , 

'  t  emp ( S i ) ' , 

'temp  (Be)  '  , 

' t  emp ( Sb ) ' , 

'temp  (As )  '  , 

'temp (As_valve )  '  , 

'is_open (Ga) ' , 

'is_open (A13) ' , 

' is_open (A14 ) ' , 

'is_open (In) ' , 

'is_open (Si) ' , 

'is_open (Be) ' , 

'is_open (Sb) ' , 

'is_open (As) ' , 

'reading (flux) ' , 

"c : / /u/ /mbe/ /recipes/ / stef an/ /log_dat a/ /out file . dat") ; 

echo("log_id  =  ",log_id); 
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/**/ 

/***************************************************************************************/ 


/*  */ 
/*  This  command  file  starts  the  pyrometer  control  program  Tsub.cmd  */ 

/*  By  launcing  it  instead  of  loading  it,  it  will  run  in  parallel  with  */ 

/*  the  rest  of  the  growth.  */ 

/*  Stefan  Svensson,  ARL  May  24  2000  */ 


/***************************************************************************************/ 

/**/ 

/**/ 

pyro_id  =  nsh ( "Tsub . cmd" ) ;  /*  Process  ID  of  the  launch  file  */ 

echo("  Launched  Tsub.cmd"); 
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■k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

k 

This  command  file  sets  a  new  substrate  temperature  over  a  time  interval  of  * 

50000  seconds  (in  other  words  it  runs  until  you  kill  it) .  * 

New  setpoints  are  entered  by  setting  the  global  variable  Tsub_new  in  an  other  * 

command  file.  Enter  the  line  load  Tsubxxx.cmd  (where  xxx  is  the  temp)  * 

Temperature  control  via  pyrometer  and  PID  routine  from  Rich  Leavitt  * 

★ 

Stefan  Svensson,  ARL  Jun  09,  2000  * 

k'k'k'kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


/**/ 

I  k  *  *  *  * 

/* 

/ 

/ 

/ 

/ 

/ 

/* 

/* 

I  k  *  *  * 

/**/ 

#include  <stdlib.h> 
#include  <signal.h> 
#include  <unistd.h> 
#include  <cells.h> 
#include  <gauges.h> 
#include  <shutters.h> 
#include  <time.h> 
#include  <mbe.h> 

/**/ 

/ *=================== 

/*=================== 

/ ^ =================== 

/**/ 


===============================================================*/ 

Declarations  and  initial  values  ==============================*/ 

===============================================================*/ 


double 

r e  c ipe_t ime 

= 

50000. 

.0;/* 

Seconds  total  time  the  cmd  file  is  running 

*/ 

double 

tOO; 

/* 

Seconds  start  time  for  cmd  file 

*/ 

double 

TP0 ; 

/* 

Pyrometer  reading  at  start  of  PID  control 

ramp 

*/ 

double 

TPstep; 

/* 

Pyrometer  step  length  in  PID  control  ramp 

*/ 

double 

TCsub; 

/* 

Thermocouple  general  substrate  temperature 

*/ 

double 

TC_test 

= 

450.0; 

:  /* 

Thermocouple  start  testing  for  pyrometer  operation 

*/ 

double 

Ramp_rate 

= 

40; 

/* 

Thermocouple  Ramp  rate  (deg/min) 

*/ 

double 

Time_step; 

/* 

Time  step  length  in  ramp 

*/ 

double 

Numb_steps 

= 

10.0; 

:  /* 

Number  of  steps  that  ramp  is  broken  into 

*/ 

int 

i; 

/* 

Counter  for  pyrometer  readings 

*/ 

int 

I_pyro 

= 

5; 

/* 

Number  of  pyrometer  readings  to  average 

*/ 

double 

newT; 

/* 

Pyrometer  New  temperature  value  used  inside 

PID  loop 

*/ 

double 

tO; 

/* 

Seconds  Start  time  for  PID  loop 

*/ 

double 

control_time; 

/* 

Seconds  duration  of  PID  loop 

*/ 

double 

Tstep 

= 

1.0; 

/* 

Seconds  wait  time  between  PID  steps 

*/ 

double 

sigma; 

/* 

PID  sigma 

*/ 

double 

delta; 

/* 

PID  delta 

*/ 

double 

Prop 

= 

72 . 0; 

/* 

PID  prop  constant 

*/ 

double 

Integ 

= 

4.7; 

/* 

PID  integration  constant 

*/ 

double 

Deriv 

= 

1.35; 

/* 

PID  derivative  constant 

*/ 

double 

TCmin 

= 

25.0; 

/* 

Thermocouple  min  allowed  value 

*/ 

double 

TCmax 

= 

850.0; 

:  /* 

Thermocouple  min  allowed  value 

*/ 

double 

TCsetp; 

/* 

Thermocouple  setpoint 

*/ 

double 

errtot 

= 

0.0; 

/* 

PID  total  error 

*/ 

double 

curr_err ; 

/* 

PID  current  error 

*/ 

double 

prev_err ; 

/* 

PID  previous  error 

*/ 

/**/ 

/*=== 


/ *=======================  check  for  setpoint  changes  ========== 

/*============================================================= 

/**/ 

tOO  =  t;  /*  Start  time  */ 

while  (  t-tOO  <  recipe_time  ) 

{ 


/* 


/**/ 

/*=== 

/*=== 


echo ("First  test  line,  Tsub_new  =  ",Tsub_new); 

if  (  Tsub_old  !=  Tsub_new  ) 

{ 


W 


Ramp  to  new  temp 
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/*======================================================================================*/ 

1**1 

echo("  Starting  pyrometer  control  to  ",Tsub_new  ); 
echo ("") ; 

TPO  =  reading (pyrometer )  ; 

TPstep  =  (Tsub_new  -  TPO ) /Numb_steps ; 

Time_step  =  fabs (TPstep) /Ramp_rate* 60 . 0 ;  /*  Time  per  step  in  sec  */ 

control_time  =  Time_step; 

i  =0; 

/*  echo("  Current  temperature  =  ",TP0);  */ 

/**/ 

while  (i  <  Numb_steps) 

{ 

i  =  i+1 ; 

newT  =  TPstep*i  +  TPO; 

echo("New  temp : newT) ; 
sigma  =  -setp (subs) *Integ*Prop/100; 
prev_err  =  reading (pyrometer ) -newT; 
set_ramp (subs,  1000)  ; 
tm  =  t; 

tO  =  t; 

while  (  t-tO  <  control_time  ) 

{ 

sleep (Tstep  -  (t  -  tm) ) ; 

curr_err  =  reading (pyrometer ) -newT; 

sigma  =  sigma+Tstep*curr_err ; 

errtot  =  errtot+Tstep*curr_err*curr_err; 

delta  =  ( curr_err-prev_err ) /Tstep; 

TCsetp  =  - (curr_err+sigma/Integ+delta*Deriv) *100/Prop; 

if  (TCsetp  <  TCmin)  TCsetp  =  TCmin; 

if  (TCsetp  >  TCmax)  TCsetp  =  TCmax; 

set_temp (subs, TCsetp) ; 

prev_err  —  curr_err; 

tm  =  tm  +  Tstep; 

tm=t ; 


/* 

/* 

/**/ 

/*== 

/*== 

/*== 

/**/ 

/* 

/**/ 


/* 


/* 


echo("  New  substrate  temp  reached  ",  newT); 
Tsub_old  =  Tsub_new; 

} 

else  */ 

{  */ 


:=============================================================*/ 

Maintain  temperature  ======================================== ^ / 

:=============================================================*/ 


echo("  Maintaining  temp  at  =", Tsub_new) ;  */ 


control_time  =5.0;  /*  Control  for  5  sec  before 

newT  =  Tsub_new;  /*  testing  for  new  Tsub 

echo ("New  temp newT) ;  */ 

sigma  =  -setp (subs) *Integ*Prop/100 ; 

prev_err  =  reading (pyrometer ) -newT; 

set_ramp (subs,  1000)  ; 
tm  =  t  ; 

tO  =  t; 


while  (  t-tO  <  control_time  ) 

{ 

sleep (Tstep  -  (t  -  tm) ) ; 

curr_err  =  reading (pyrometer ) -newT; 

echo ("newT  =  ",newT);  */ 
sigma  =  sigma+Tstep*curr_err; 

errtot  =  errtot+Tstep*curr_err*curr_err; 

delta  =  ( curr_err-prev_err ) /Tstep; 

TCsetp  =  - (curr_err+sigma/Integ+delta*Deriv) *100/Prop; 


if  (TCsetp  <  TCmin)  TCsetp  =  TCmin; 


*/ 

*/ 
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Tsub.cmd 


if  (TCsetp  >  TCmax)  TCsetp  =  TCmax; 
set_temp (subs, TCsetp)  ; 
prev_err  =  curr_err; 
tm  =  tm  +  Tstep; 

tm  =  t; 

} 
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/**/ 

/******** 

/* 

/*  This 
/* 

/******** 

/**/ 


****************************************************************************** I 

*/ 

command  file  sets  the  substrate  temperature  during  growth  */ 

Stefan  Svensson,  ARL  Apr  21,2000  */ 

****************************************************************************** i 


Tsub_new  =  550; 

tm  =  t  ; 

echo ("Tsub_new=  ",  Tsub_new) ; 


19 


Appendix 


Kill_pyr.cmd 


/**/ 

/***************************************************************************************/ 


/*  */ 

/*  This  command  file  kills  the  pyrometer  control  program  Tsub.cmd  */ 

/*  */ 

/*  Stefan  Svensson,  ARL  May  24  2000  */ 


/***************************************************************************************/ 

/**/ 

/**/ 

#include  <signal.h> 

echo ("pyro_id=  ",pyro_id); 

kill  (pyro_id, SIGKILL)  ; 

echo("  Pyrometer  control  finished."); 
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/************************************************************************/ 


/*  */ 

/*  This  command  procedure  kills  the  data  log  */ 

/*  The  data  is  stored  in  the  file  */ 

/*  C : /u/mbe/recipes/stef an/log_data/outf ile . dat  */ 

/*  */ 


/************************************************************************/ 

#include  <signal.h> 
echo ("sigkill=", SIGKILL) ; 
echo ("log_id=",  log_id)  ; 
kill (log_id, SIGKILL)  ; 

echo("  Stdlogl  finished.  Data  is  in  c : \\u\\mbe\\recipes\\stefan\\log_data\\outf ile . dat") ; 
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